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摘 要 :本文 从 SMARTS 数据 库 和 之 前 文献 中 , 收集 了 次 变 体 36 454. 3(2 = 0. 859) 光学 波段 


B、V、R 和 红外 波段 J、 KK 的 数据 , 时间 从 2008 年 6 月 份 到 2017 年 7 月份 之 间 近 十 年 . 分别 
获得 了 光学 B 波 段 891 组 、V 波 段 855 组 、R 波 段 877 组 、J 波段 860 组 和 K 波 段 751 组 共 
4234 组 数据 . 使 用 这 些 观测 数据 , 用 LombScargle Periodogram(LSP) 方 法 和 加 权 小 波 Z 变换 
法 (WWZ) 研究 了 各 个 波段 的 周期 , 结果 表明 : (1)3C 454.3 在 光学 和 红外 波段 都 有 454 天 的 周 
期 光 变 ; (2) 五 个 波段 存在 454 天 的 长 周期 光 变 ,由 此 得 到 30 454.3 FS RAAE 


为 : Nf ~% 0.24 x 10°N 和 辐射 区 半径 为 : R = 3.68 x 10°Xm; (3) 通 过 长 周期 光 变 分 析 ， 


FR AH EK 30 454. 3 在 2021 年 6 月 左右 ,将 再 次 爆发 。 
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耀 变 体 的 多 波段 的 光 变 时 标 是 一 个 重要 的 物理 参数 , 焰 变 体 的 长 周期 光 变 时 标 可 以 帮助 我 们 
研究 天 体 的 中 心 黑 洞 质量 , 内 部 结构 , 辐射 区 域 等 问题 , 目前 研究 此 类 问题 的 方法 很 多 , 有 结 
构 函 数 法 , 离散 相关 函数 法 (DCF) , period4 方 法 ， 功 率 谱 (PSD), Jurkevich 方法 , 加 权 小 波 Z 


[] » op [2] Wy y 
变换 法 (WWZ) ,LombScargle Periodogram(LSP) 方 法 .由 于 天 体 的 光学 波段 和 红外 波段 的 


数据 观测 受 天 气 设备 等 多 种 因素 的 影响 , 观测 数据 在 时 间 序列 上 不 连续 , 在 光 变 周 期 的 研究 
中 引入 了 较 大 的 误差, 结果 的 精确 度 一 直 是 该 研究 领域 的 热点 问题. 在 文献 [5] 中 对 这 个 天 体 
的 光学 和 红外 波段 的 数据 都 采用 了 功率 谱 (PSD) 分 析 方 法 , 也 得 到 了 与 本 文 相同 的 结果 , 但 是 
在 功率 谱 方法 中 要 求 数据 为 等 间隔 , 所 以 在 文献 [5] 中 对 相 邻 的 数据 进行 了 平均 插值 处 理 , 这 
样 可 能 会 对 原始 数据 的 真实 性 有 所 损害 , 还 引入 了 人 为 主观 的 因素 , 而 且 还 可 能 会 带 来 其 他 


n ZHANG PR. AICASMARTSHI CHR [5] HERT ekc 454.3 Z = 0.859) © 


的 准 同时 性 光学 波段 近 十 年 的 数据 .采用 了 LSP 方 法 来 对 耀 变 体 3C 454. 3 的 非 均 匀 数 据 进 行 
分 析 , 在 LSP 方 法 中 则 不 需要 对 数据 进行 插值 处 理 , 但 是 在 LSP 方 法 的 功率 谱 结果 中 可 能 会 出 
现 虚 假 的 峰值 , 为 了 进一步 的 评估 功率 谱 结 果 峰 值 的 可 信和 度 , 通过 模拟 大 量 光 变 曲线 , 并 采用 
LSP 方 法 计算 模拟 光 变 曲线 的 功率 谱 , 对 模拟 光 变 曲线 的 功率 谱 结果 进行 抽样 , 估计 出 相应 波 
段 的 置信 度 , 通过 置信 度 可 以 评估 所 出 现 峰 值 的 可 靠 性 , 我 们 就 可 以 得 出 更 为 准确 的 准 周期 ， 
并 且 在 分 析 结 果 中 我 们 看 到 , 所 得 到 的 周期 就 只 一 个 , 而 没有 其 他 的 周期 干扰 , 将 这 种 模拟 大 
量 光 变 曲线 的 方法 称 作 为 蒙特 卡 洛 模拟 方法 , 为 了 验证 蒙特 卡 洛 模拟 方法 所 得 到 结果 的 可 靠 
性 和 可 行 性 , 我 们 采用 这 五 个 波段 的 数据 , 用 加 权 小 波 Z 变 换 WWZ) 方 法 进行 长 周期 光 变 分 析 ， 


也 得 到 相同 的 结果 , 并 通过 张 雄 等 人 文献 中 ”标准 差 公式 , 计算 相应 波段 的 置信 和 度 , 得 出 了 可 


靠 的 准 周期 , 从 而 验证 了 蒙特 卡 洛 模拟 方法 所 得 到 结果 的 可 行 性 和 准确 性 . 在 本 文中 首次 将 
耀 变 体 3C 454. 3 光学 B, V, R 和 红外 J, kK 波 段 的 离散 数据 应 用 加 权 小 波 Z 变 换 法 (WWZ), 和 蒙特 卡 
洛 模拟 方法 来 研究 相应 波段 的 光 变 周期 , 为 泡 变 体 的 长 周期 光 变 的 研究 找到 了 一 种 新 方法 ， 
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通过 WWZ 和 蒙特 卡 洛 模拟 方法 可 以 很 好 的 得 到 非 平稳 信号 的 准 周期 , 大 大 提高 了 相应 周期 的 
计算 精度 , 从 而 得 到 更 加 准确 的 准 周期 . 在 这 两 种 方法 中 我 们 得 到 了 五 个 波段 都 存在 约 为 454 


天 的 准 光 变 周期 . 从 长 周期 的 下 
而 为 活动 星系 核 (AGN) 的 物理 


米 


完 中 我 们 可 以 得 到 3C 454. 3 中 心 黑洞 质量 和 辐射 区 半径 , 从 
模型 的 研究 提供 重要 参数 . 
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观测 数据 和 研究 : 


如 图 1 所 示 的 为 来 自 于 SMARTS 数 据 库 和 文献 [5] 的 耀 变 体 3C 454. 3 光学 B, R, V 波 段 和 红外 


J,K 波 段 的 数据 所 给 出 的 光 变 曲线 , 横 轴 为 儒 略 日 (MJD), 纵 轴 为 星 等 值 (MAG) : 


54500 55000 55500 56000 


曲线 , 紫色 的 为 红外 J 波段 的 光 变 曲线 , 呈 


56500 57000 57500 58000 


MJD (day) 
All: 红色 的 为 光学 B 波 段 的 光 变 曲线 , 检 色 的 为 光学 V 波 段 的 光 变 曲线 , 蓝 色 的 为 光学 R 波 段 的 光 变 


色 的 为 红外 kK 波段 的 光 变 曲线 


The red, orange, and blue, line represents a light curve in the optical B-Band, V-Band, 


R-Band, the pupre, and black line represents a light curve in the infrared J-Band, and 


K-Band, respectively. 


采用 流量 与 星 等 的 转化 关系 , 将 以 上 的 五 个 波段 的 星 等 值 转化 为 流量 值 , 公式 如 下 : 


F =10 25 x3640Jy (1) 
FE 为 每 一 个 波段 所 对 应 的 流量 值 ， 单 位 为 Jy， 
1.1 光 变 曲线 的 长 周期 分 析 : 


从 SMARTS 数据 库 和 文献 [5 中 获得 星 等 的 


m 为 所 对 应 波段 的 星 等 值 . 


数据 之 后 , 通过 流量 与 星 等 的 转化 公式 得 到 了 


[6] 


各 个 波段 的 流量 数据 . 接 下 来 本 文 将 采用 蒙特 卡 洛 模拟 “和 加 权 小 波 2 变换 法 (WWZ) 来 计算 各 
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个 波段 的 周期 和 各 自 的 置信 度 . 即 用 如 图 1 所 示 的 光 变 曲线 数据 获得 流量 数据 之 后 , 采用 LSP 
方法 得 到 相应 的 光 变 周期 和 蒙特 卡 洛 模 方 法 得 到 相应 的 置信 度 估 计 , 为 了 检验 蒙特 卡 洛 模拟 
方法 置信 度 评估 的 可 行 性 和 准确 性 , 我 们 通过 了 加 权 小 波 Z 变 换 法 (WWZ) , 计算 了 相对 应 波段 
的 光 变 周期 , 通过 标准 差 公式 , 我 们 计算 了 相应 的 加 权 小 波 Z 变 换 周期 的 置信 和 度 . 


1.2 应 用 LombScargle Periodogram(LSP) 方 法 分 析 长 周期 光 变 : 

计算 隐藏 在 噪声 中 的 周期 性 信号 是 天 文 数据 时 间 序 列 分 析 中 的 一 个 重要 目标 .LSP 方 
法 可 以 对 非 均 匀 采 样 的 时 间 序 列 周期 图 进行 相位 修正 处 理 , 能 够 在 一 定 范围 内 对 非 均匀 
采样 的 时 间 间 隔 引 起 的 误 报 周期 进行 修正 , 并 且 LSP 方 法 不 要 求 数据 是 均匀 的 时 间 序 列 ， 
不 用 对 其 进行 插值 , 从 而 减少 了 人 为 因素 对 数据 的 干扰 . 因此 , LSP 方法 能 很 好 的 寻找 隐 
藏 在 噪声 中 的 准 周 期 振荡 光 变 . 现在 假定 有 这 么 一 组 时 间 序 列 x(t), (7 = 1, 2, 


ee A), SUS TaD AP AES LSP 的 功率 谱 由 以 下 的 公式 所 给 出 , Px Co, ) 以 角 


Z (@, = 2zmy)) 为 变量 的 功率 谱 函数 , 基本 形式 如 下 : 


SR Flos o,(¢, —r) Yb) —Flsin@,(¢,—7) 
P(w,)=~ E + &§ m O) 


2 N N 5 
> cos’ w(t, -7) > sin @,(t, —7) 
j=l j=l 


N 
这 里 的 = 二》 al 为 时 间 序列 的 平均 值 , Y 为 数据 的 个 数 ，r 为 相应 的 时 间 的 相位 修 
j=l 


iE: 


tan(2@,T) = (3) 


对 于 图 1 的 光 变 曲线 是 一 个 非 均 匀 的 时 间 序列 , 要 得 到 里 面 所 隐藏 真实 周期 是 非常 困难 
的 .采用 以 上 的 LSP 等 式 分 析 真 实数 据 的 功率 谱 结果 , 但 是 由 于 这 些 结果 中 会 出 现 虚假 
峰值 , 使 得 我 们 很 难 判别 真实 的 周期 频率 , 所 以 为 了 解决 这 个 问题 我 们 需要 对 这 些 功率 
谱 的 峰值 进行 置信 度 的 评估 


1. 3 蒙特 卡 洛 模拟 方法 和 置信 度 评估 ; 
AGN 的 功率 谱 一 般 呈 红 噪 声 守 律 分 布 名 即 : P ec £7? 为 了 获得 功率 谱 指 数 ui ， 对 所 得 


到 的 功率 谱 结果 取 对 数 , 并 进行 一 元 线性 回归 拟 合 , 可 获得 功率 谱 指数 q 值 .图 2 中 的 横 
坐标 频率 f 的 对 数值 , 纵 坐 标 为 功率 谱 结果 的 对 数值 , 红色 直线 为 相应 的 一 元 线性 拟 合 
直线 , 图 中 直线 的 斜率 为 相应 的 功率 谱 指数 q . 
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图 2 光学 B, R, V 波段 和 红外 J, K 波段 的 功率 谱 对 数 坐标 拟 合 


Fig. 2 Logarithmic coordinate fitting of power spectrum in optical B,R,V bands and infrared J,K bands 
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通过 以 上 的 一 元 线性 回归 拟 合 , RAET NAAA R KARE a. 基于 功 
率 谱 指数 a 和 真实 数据 , 模拟 了 5000 条 光 变 曲线 , 对 模拟 的 光 变 曲线 使 用 LSP 方法 逐一 进 
行 计算 , 获得 模拟 光 变 曲线 的 功率 谱 结果 , 并 对 这 些 模拟 光 变 曲线 的 功率 谱 结果 进行 置信 
度 抽样 , 可 得 到 各 个 波段 不 同 的 置信 度 曲 线 , 最 终 将 真实 数据 的 功率 谱 和 置信 度 曲 线 整合 ， 
即 可 得 到 准确 的 周期 频率 . 


1.3 LSP 方法 的 置信 和 度 评估 结果 : 
3 中 分 别 为 光学 B, R, V 波段 , 和 红外 J, K 波段 的 蒙特 卡 洛 模 拟 LSP 方法 结果 , 图 3 中 纵 
轴 为 相应 波段 的 功率 谱 值 , 横 轴 为 相应 的 频率 了 ， 单 位 为 d-, 红色 曲线 代表 了 真实 数据 
的 功率 谱 结果 , 黑色 的 曲线 代表 了 置信 度 为 : 99.7%, 紫色 的 曲线 代表 了 置信 度 为 : 
99% , EER HEARSE 7 BUREN: 95%, 我 们 取 超 过 置信 度 曲线 99.79% 的 峰值 为 所 对 
应 的 周期 频率 , 如 图 中 黑色 箭头 所 指 的 位 置 . 

图 3 光学 B,RV 波段 和 红外 JK 波段 的 置信 度 评估 结果 


Fig. 3 Confidence assessment results of optical B,R,V bands and infrared J,K bands 
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由 图 3 我 们 看 到 了 在 光学 B,R,V 和 红外 J,K 波 段 都 只 存在 一 个 很 明显 的 峰值 超过 了 
99.7% ,如 图 3 中 箭头 所 指 的 位 置 , 相应 的 峰值 频率 保留 四 位 小 数 分 别 都 为 


F = 0. 0022(d 7). 所 以 获得 光学 和 红外 波段 都 存在 约 为 454 天 的 周期 ,通过 这 样 的 方 
法 我 们 得 到 了 与 之 前 文献 [5] 中 相同 的 结果 , 并 且 通 过 这 种 方法 不 会 引入 其 他 干扰 结果 . 
为 了 检验 以 上 方法 的 可 行 性 和 准确 性 , 我 们 通过 了 加 权 小 波 2 变换 法 (WWZ) ,计算 了 相 
应 波段 的 周期 . 由 于 标准 差 和 置信 度 的 关系 , 我 们 计算 了 相应 的 加 权 小 波 Z 变换 功率 谱 


n [4] 


的 置信 度 , 本 文中 采用 通过 之 前 文献 “中 所 给 出 的 正太 分 布 的 标准 差 公式 . 


2 加 权 小 波 Z 变换 法 (WWZ) 分 析 长 周期 光 变 : 

在 经 典 的 时 间 频 率 分 析 中 一 般 采 用 傅 里 叶 变 换 和 小 波 分 析 方 法 , 但 是 由 于 天 文 数据 往往 
都 是 非 等 间距 的 离散 数据 , 在 使 用 傅 里 叶 变 换 和 小 波 变 换 来 分 析 天 文 数据 时 需要 对 数据 进行 
i 值 处 理 变 成 均匀 的 数据 , 这 样 对 数据 的 真实 性 由 很 大 的 损害 , 并 且 在 使 用 传 里 叶 变换 来 处 
这 些 非 等 间距 的 数据 时 可 能 还 会 出 现 伪 周 期 , 所 以 在 基于 小 波 变换 的 基础 上 Foster“ 提 
了 加 权 小 波 Z 变换 (WWZ), 这 样 不 仅 可 以 更 加 有 效 的 处 理 非 等 间距 的 时 间 序 列 的 周期 , 还 可 
以 一 定 程度 上 的 反应 出 周期 的 稳定 性 . 加 权 小 波 2 变换 是 将 时 间 序 列 投影 到 三 个 正 交 归 一 的 


TERME, AC) = 1, gE) = coslat — 1,)1,¢@,) = sinlot - n) 1, ERZE 


oar 


FE [ii 


做 了 统计 加 权 w，= expl-co(t, 一 zt)], 将 数据 的 不 均匀 通过 权重 调节 , 避免 了 周期 分 
析 时 受到 数据 过 密 的 影响 , 其 中 使 用 的 母 函 数 为 Morlet 小 波 ,WWZ 的 定义 为 : 


(Neff 一 3)Vy 


WWZ= 
2(Vz — Vy) (4) 


且 这 个 等 式 的 分 子 分 母 分 别 满足 自由 度 为 W_ ，- 3 和 2 的 F 分布. 其 中 /Vj 为 有 效 数 


据点 个 数 , ,分别 为 观测 数据 和 模拟 函数 的 加 权 变 量 , 即 : 


bo expl-2co (t, zi 
a >) exp [-2cw*(t, -r,)*] 


N 
(5) 
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4 a 4 (6) 
Yast) |> oy t) | 
V i 

” 2 o, 2 o, 


a a | (7) 
Hc 为 衰减 因子 ,7 为 时 移 , @ 为 尺度 因子 , @, (A =1,2,3......... n) 为 对 应 的 测试 频 


率 
PAS 


图 4 光学 B,R,V 波段 和 红外 ],K 波段 的 WZ 分 析 结果 和 相应 的 周期 频率 . 


Fig. 4 the WWZ analysis results and corresponding to periodic frequency in the 
optical B,R, V bands and infrared J and K bands. 
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从 以 上 的 WZ 结果 图 中 可 以 得 到 在 光学 B, R, V 和 红外 J,K 波 段 左边 彩 图 红色 部 位 对 应 了 
右 图 的 功率 谱 曲 线 所 出 现 的 峰值 频率 的 位 置 为 相应 的 周期 频率 的 出 现 的 位 置 , 通过 白 骂 


j 
声 检验 采用 “， be Ji- r) , 这 个 等 式 对 于 正 态 分 布 的 总 体 是 成 立 的 , 等 


AH ô, = 1 一 克 为 相应 波段 的 每 一 个 功率 值 , 17 相应 波段 功率 的 平均 值 , M 为 相应 的 离 


散 功 率 谱 个 数 , 通过 这 个 等 式 计算 以 上 各 个 波段 的 功率 谱 的 标准 差 , 由 于 正太 分 布置 信 
度 与 标准 差 的 关系 , 我 们 可 以 获得 各 个 波段 的 功率 谱 的 置信 度 , 采用 3 售 标准 差 所 对 应 
的 置信 度 99. 7%, 所 以 图 4 中 右 图 的 蓝 色 直 线 取 3 倍 的 标准 差 值 . 通过 以 上 的 结果 各 个 波 
段 周期 频率 出 现 的 位 置 都 为 = 0. 0022A), 所 以 光学 和 红外 都 存在 约 为 454 天 的 


周期 . 并 且 其 置信 和 度 都 超过 了 99.7%, 这 个 结果 与 上 面 的 蒙特 卡 洛 模拟 LSP 方法 的 结果 
一 致 . 综 上 我 们 得 到 了 炮 变 体 3C 454. 3 在 光学 B,R,V 和 红外 J,K 波 段 的 周期 都 为 454 R, 


大 约 为 1. 244 年 . 
3. 3C 454.3 中 心 黑洞 质量 和 辐射 区 域 半 径 的 估计 : 
3.1 中 心 黑 洞 质量 估计 : 
中 心 黑 洞 质量 的 研究 也 是 人 们 极为 感 兴趣 的 , 汾 变 体 中 心 有 一 个 超大 质量 黑洞 ,很 多 辐射 现 


象 是 由 中 心 黑洞 所 引起 的 . 通过 长 时 标 光 变 能 够 获得 中 心 黑洞 的 质量 ,假设 由 薄 吸 积 盘 理 论 


[7] 


引起 的 长 周期 光 变 , 那么 由 之 前 文献 中 “的 公式 : 


fy, = 4.5202 ME? yrs (8) 


等 式 中 ，p 为 HS eM, = MOM, 为 中 心 黑洞 质量 约 化 单位 , 吸 积 率 M 和 


8 
~ MM AK ES Bu, Mu = 0.5， M =M, |e Ms 是 爱 丁 顿 吸 积 率 ，& 是 吸 


RACK. 爆发 时 间 间 隔 主 要 取决 于 .Horiuchill 和 Kato(1990) 建 议 到 4 = 0. 5 时 ,磁场 


© 的 逃逸 速率 会 更 低 . 此 时 , 热 极限 周期 时 间 为 上 = 2 


burts “ 


© tye = 9.0855 Mg” yrs (9) 


west B= B/0.1, PUMP = 0.1, u = 0.5, M=0.2Mc, 时 我 们 采用 五 个 波 


段 的 1. 244 年 的 周期 得 到 中 心 黑洞 的 质量 约 为 : = 0.24 x 1L0°U, . 这 个 结果 是 由 于 利 


用 薄 吸 积 盘 理 论 的 分 析 方 法 , 没有 考虑 到 黑洞 的 自 旋 的 影响 , 所 以 比 一 般 的 平 谱 射电 类 星体 
的 中 心 黑洞 质量 偏 小 . 

r 3.2 辐射 区 半径 : 

© 对 于 兆 变 体 的 长 周期 光 变 的 现象 , 目前 还 没有 最 佳 的 物理 模型 解释 , 现 有 的 理论 为 双 黑 洞 


-He [8] » st Hs ate He rey 5 = a ` ` 
模型 ”， 薄 盘 的 热 不 稳定 性 模型 ”, 和 螺旋 喷 流 等 模型 ”. 如 果 是 薄 盘 不 稳定 性 造成 的 , 可 


从 所 搜集 的 数据 中 获得 各 个 波段 的 星 等 变化 情况 : 
表格 1. 每 个 波段 的 幅 值 变化 


Table 1. the variation of magnitude in each band. 


Wavebands MaximUm Minimum Difference 
B 17.337 13.946 3.373 
V 16.756 13.448 3.308 
R 16.585 12.97 3.615 
J 15.092 11.012 4.08 
K 15.017 9.188 5.829 


由 以 上 表格 1, 可 以 看 见 ,在 光学 和 红外 波段 都 具有 较为 剧烈 的 变化 , 并 且 红外 波段 的 变化 
幅度 比 光学 波段 的 变化 更 加 剧烈 . 采用 薄 盘 不 稳定 性 分 析 长 周期 “， 通 过 长 周期 “可 


以 确定 出 热 不 稳定 性 产生 的 区 域 , 等 式 如 下 : 


1 
T =3.2x10* B’ (x?(x-1)M,) yrs (10) 


ERP x = R/R, R, 为 相应 的 施 瓦 西 半径 , BON Bi RL, M = M/10'Mo. 这 里 我 们 


取 :7 = 1.244yrs, 中 心 天 体 黑洞 质量 为 : M = 0.24 x 10° UL, . 所 以 我 们 得 到 辐射 区 
域 半 径 为 : R = 3.68 x 10? fom. 


4 . 结果 


YEAR BS 3C 454. 3 的 光学 波段 和 红外 波段 的 研究 中 , 本 文采 用 了 LSP 方法 得 到 了 各 个 波段 

的 周期 并 评估 了 所 出 现 的 峰值 频率 的 置信 度 , 得 到 了 更 为 可 靠 的 结果 , 并 采用 了 WZ 方法 

验证 蒙特 卡 洛 模拟 评估 置信 度 方 法 的 结果 , 我 们 发 现 了 : 

(1) 两 种 方法 都 在 光学 B,R,V 波段 和 红外 于 波段 发 现 了 454 R, 约 为 1. 244 年 的 光 变 周 
期 , 通过 蒙特 卡 洛 模拟 评估 置信 度 方 法 , 我 们 清楚 的 看 到 了 对 于 五 个 波段 都 只 有 一 个 

峰值 频率 超过 了 99. 7%, 并 且 没 有 引入 其 他 干扰 结果 . 而 为 了 进一步 的 验证 这 种 方法 的 

] 行 性 , 我 们 对 这 五 个 波段 采用 了 WZ 的 分 析 方 法 , 也 得 到 了 与 蒙特 卡 洛 模拟 方法 一 

致 的 结果 , 都 为 454 天 , 约 1. 244 年 . 

> (2) 对 于 中 心 的 超大 质量 黑洞 , 通过 之 前 文献 所 给 出 的 长 周期 与 黑洞 质量 的 公式 , 结合 

~ 本 文中 的 五 个 波段 共同 的 长 周期 数据 我 们 得 到 了 中 心 黑洞 的 质量 


为 : M x 0. 24 x LOM, ,这 个 结果 由 于 是 利用 薄 吸 积 盘 理论 的 分 析 方 法 , 没有 考 上 
到 黑洞 的 自 旋 的 影响 , 所 以 比 一 般 的 平 谱 射电 类 星体 的 中 心 黑 洞 质量 偏 小 , 辐射 区 半 
H: R ~ 3. 68 x 10 Kan ,计算 所 得 到 的 辐射 区 半径 相 比 之 前 文献 ”所 给 出 的 和 


z| 


Er 


Hal 


É 射 区 半径 在 相应 的 施 瓦 西 半 径 下 要 小 , 由 于 之 前 文献 no 中 所 采用 的 黑洞 质量 要 比 本 


一 文中 所 用 的 黑洞 质量 大 , 所 以 此 处 我 们 所 得 的 辐射 区 半径 是 合理 的 . 

© (3) 通过 以 上 的 计算 , 我 们 预期 兆 变 体 3C 454. 3 在 2021 年 6 AAA, 将 再 次 爆发 . 我 们 将 

用 云南 天 文 台 丽 江 的 2.4 米 光学 望远镜 做 进一步 实测 验证 我 们 的 结果 . 为 用 蒙特 卡 洛 
模拟 方法 研究 次 变 体 的 长 周期 光 变 寻 求 观测 证 据 . 

致谢 :真诚 致谢 SMARTS 团队 和 吴 月 承 等 人 在 文献 [5 中 所 提供 的 光学 和 红外 的 星 等 数据 . 
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Studying the properties of long-period light variation about blazar 3C 454.3 with 

the methods for Lomb-Scargle Periodogram 

Lu Lin, Zhang Hao-Jing Zhang Huan, Ma kai-Xuan 

(Department of Physics, Yunnan Normal University, Kunming 650092, China) 
Abstract: In order to study some feature of optical and infrared radiation in this paper. The data 
are collected from SMARTS dataset and other literature about blazar 3C 454.3(z=0.589), which 
are the optical bands B, V, R and infrared bands J, K. The data covered from June 2008 to July 
2017 about decade, were obtained with 891 groups in B-band, 855 groups in V-band,877 groups 
in R-band, 860 groups in J-band, and 751 groups in K-band, there are 4,234 points of data. Using 
the observed data, we employ LombScargle Periodogram (LSP) and The Weighted Wavelet 
Z-Transform(WWZ) method to study the possible periods in each band, the results show that: (1) 
There is a period of 454 days in both optical and infrared bands of 3C 454.3 ; (2) Using a period of 
454 days, we obtain the mass of black hole about the 3C 454.3, approximate 


to:N ~ 0.24 x 10°M, and the radius of radiation region: R ~ 3.68 x 10° Am (3) Through 


these analyses, the blazar 3C 454.3 will burst out again around June 2021, 


Key words: Monte Carlo Simulation, The Weighted Wavelet Z-Transform, 
the mass of black hole, Radius of radiation region 


